INSTABILITY OF TOROIDAL MAGNETIC FIELD 
IN JETS AND PLERIONS 

Mitchell C. Begelman 1,2 

JILA, University of Colorado and National Institute of Standards and Technology, 

Boulder, CO 80309-0440 

Submitted to Ap. J., April 28, 1997; in revised form, July 16, 1997; accepted August 2 1997. 

ABSTRACT 

Astrophysical jets and pulsar-fed supernova remnants (plerions) are expected to de- 
velop highly organized magnetic structures dominated by concentric loops of toroidal field, 
Bff,. It has been argued that such structures could explain the polarization properties of 
some jets, and contribute to their lateral confinement through magnetic tension forces. A 
concentric toroidal field geometry is also central to the Rees-Gunn model for the Crab 
Nebula — the archetypal plerion — and leads to the deduction that the Crab pulsar's 
wind must have a weak magnetic field. Yet this kind of equilibrium between magnetic and 
gas pressure forces, the "equilibrium Z-pinch" of the controlled fusion literature, is well 
known to be susceptible to disruptive localized instabilities, even when the magnetic field is 
weak and/or boundary conditions (e.g., a dense external medium) slow or suppress global 
modes. Thus, the magnetic field structures imputed to the interiors of jets and plerions 
are unlikely to persist for very long. 

To determine the growth rates of Z-pinch instabilities under astrophysical conditions, 
I derive a dispersion relation that is valid for the relativistic fluids of which jets and ple- 
rions may be composed, in the ideal magnetohydro dynamics (MHD) limit. The dominant 
instabilities are kink (m = 1) and pinch (m = 0) modes. The former generally dominate, 
destroying the concentric field structure and probably driving the system toward a more 
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chaotic state in which the mean field strength is independent of radius (and in which re- 
sistive dissipation of the field may be enhanced). I estimate the timescales over which the 
field structure is likely to be rearranged and relate these to distances along relativistic jets 
and radii from the central pulsar in a plerion. 

I conclude that the central tenet of the Rees-Gunn model for the Crab Nebula — the 
existence of a concentric toroidal field well outside the pulsar wind's termination shock 
- is physically unrealistic. With this assumption gone, there is no dynamical reason to 
conclude that the magnetic energy flux carried by the pulsar wind is much weaker than the 
kinetic energy flux. Abandoning the principal conclusion of Rees & Gunn would resolve a 
long-standing puzzle in pulsar wind theory. 

1. INTRODUCTION 

The poloidal (i.e., radial: B z ) and toroidal (B^) components of magnetic field frozen 
into a perfectly conducting outflow behave differently as the flow expands. If the fluid in a 
jet with cylindrical radius r(z) <C z and speed v(z) expands laterally as it moves outward 
along the z-axis, the poloidal field varies as r~ 2 while the toroidal field varies as r~ l v~ x 
(Begelman, Blandford, & Rees 1984). Thus the toroidal field should eventually become 
dominant. The toroidal component of the field is expected to be even more dominant in 
the wind from a pulsar, since it is amplified initially by the pulsar's spin. In the wind 
from the Crab pulsar, for example, the toroidal and poloidal field components should be 
comparable near the light cylinder radius, 2 x 10 8 cm. By the time the wind reaches the 
reverse shock (and inner boundary of the visible nebula), at about 10 17 cm, B^ is expected 
to be nearly 9 orders of magnitude stronger than B z . 

Theorists have made much of the possible astrophysical importance of organized field 
structures dominated by concentric loops of toroidal magnetic field. First, such field struc- 
tures can have important dynamical consequences, for example, providing substantial pres- 
sure confinement in a jet (Benford 1978; Chan & Henriksen 1980; Bicknell & Henriksen 
1980; Eichler 1993), via magnetic tension. Second, the presumed existence of a highly or- 
ganized field has been used to place tight constraints on otherwise unobservable properties 
of the Crab pulsar's wind. By assuming a concentric toroidal field in the Crab Nebula, 
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Rees & Gunn (1974) were able to deduce that the kinetic energy flux in the wind must be 
2-3 orders of magnitude larger than the Poynting flux, at the radius at which the wind 
encounters the reverse shock. This empirical conclusion has been confirmed and refined 
by later analyses (Kennel & Coroniti 1984a; Emmering & Chevalier 1987; Begelman & 
Li 1992) but has proven hard to reconcile with theoretical models of pulsar winds, which 
predict comparable fluxes of kinetic and electromagnetic energy (Michel 1969; Kennel, Fu- 
jimura, & Okamoto 1983; Coroniti 1990; Arons 1992; Begelman & Li 1994; Tomimatsu 
1994; Michel 1994; Melatos & Melrose 1996). 

The applications cited above depend critically on the assumption that the field lines 
remain concentric within the Crab Nebula, or within an expanding jet. But there are good 
reasons to question the plausibility of this assumption. Plasma physicists have known for 
more than 40 years that simple field configurations consisting of concentric circular loops 
of field are generically unstable, even when the magnetic pressure is weak compared to the 
gas pressure (Kruskal & Schwarzschild 1954; Kadomtsev 1966). Such systems, called "Z- 
pinches" in the plasma physics literature (because the current flows in the z— direction), 
are subject to short- wavelength modes that cannot be suppressed by a judicious choice of 
boundary conditions, such as a very dense external medium or a highly supersonic flow 
velocity (e.g., Konigl & Choudhuri 1985). Even if the poloidal field were amplified by 
shear, as seems to occur in some jets, the poloidal field structure would lack the large- 
scale coherence — in contrast to — necessary to suppress kink (m = 1) and pinch 
(m = 0) instabilities through tension forces. In particular, modes with an m = 1 angular 
dependence (i.e., oc e 1 ^) will destroy the concentric symmetry, leading to a much less 
regular field structure. 

In this paper, I discuss the effects of short-wavelength Z-pinch instabilities as they 
apply to astrophysical jets and pulsar-fed nebulae (plerions) such as the Crab Nebula. By 
"short" wavelength, I mean wavelengths that are much shorter than any of the length 
scales associated with the unperturbed medium, such as the width of the jet or radius 
of the pier ion. These modes are driven by forces due to the curvature of field lines and 
magnetic and gas pressure gradients in the background equilibrium. They are distinct from 
the Kelvin- Helmholtz modes, driven by shear in the velocity field, which have been studied 
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extensively in connection with jets. These local modes are also distinct from the global 
pinch and kink modes that can affect magnetized jets under certain conditions (Chan & 
Henriksen 1980; Eichler 1993). I estimate the growth rates of the most unstable Z-pinch 
modes, and assess how these modes might reorganize the structure of the magnetic field, 
given the available evolutionary timescales. 

I conclude that, while the magnetic field in the Crab Nebula and many jets may 
remain aligned with the toroidal plane, the large-scale concentric structure of the field is 
probably destroyed by instabilities. This conclusion, if true, would obviate Rees & Gunn's 
dynamical argument that the kinetic energy flux must exceed the Poynting flux in the 
Crab pulsar wind. 

The plan of this paper is as follows. Since the properties of Z-pinch instabilities seem 
not to be very familiar to the astrophysical community, I present heuristic derivations of 
the m = and m = 1 stability criteria in § 2. Readers may also wish to consult standard 
textbooks (especially Bateman 1978; also Krall & Trivelpiece 1973; Chen 1974; Jackson 
1975) and monographs (e.g., Freidberg 1982) for further background on these well-studied 
instabilities, particularly in the context of controlled fusion devices. The full local (i.e., 
short- wavelength) dispersion relation is derived in § 3, taking special relativistic effects into 
account. I am not aware of other fully relativistic calculations of the Z-pinch dispersion 
relation in the literature. (In deriving the dispersion relation, I include perturbations of 
a small, uniform B z . This can either enhance or detract from stability, depending on the 
helicity of the mode.) Important consequences of the instability, notably growth rates and 
nonlinear development, are discussed in § 4. Section 5 deals with applications to jets and 
the Crab Nebula. I summarize the main results in § 6. 

2. HEURISTIC DERIVATION OF STABILITY THRESHOLD 

Consider a cylindrically symmetric MHD equilibrium with a magnetic field of the form 

B = B(r)0 + bz, (2.1) 

where b is a constant <C B(r). We assume that b is independent of r so that it does not 
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figure in the equilibrium condition, 

where p is the gas pressure. Despite its irrelevance to the equilibrium state, the inclusion 
of a small B z component has an important effect on stability. It can either suppress or 
enhance instability at sufficiently short wavelengths, and dramatically affects the formal 
stability threshold for m = (axisymmetric) modes, as we will show below. Strictly 
speaking, a field configuration with the form (2.1) is a special case of a "screw pinch" 
(Freidberg 1982), because individual field lines are cylindrical helices; however, we will 
continue to use the term Z-pinch to express its principal characteristic. 

To provide some physical insight into a Z-pinch, we follow what happens when a 
narrow ring of magnetic flux is displaced slightly against an unperturbed background. For 
ease of visualization, we restrict our analysis to m = 1 (kink) and m = (pinch) modes. 

Consider first an m = 1 displacement, which corresponds to moving a flux loop side- 
ways in the r — <p plane (neglecting the helical rise) without changing its radius or circular 
shape. If we make this change slowly enough (simulating conditions near marginal stabil- 
ity), then the gas pressure will remain uniform around the loop; since the volume of the 
loop is unchanged to first order, the gas pressure is also unchanged, 

pi=Po- (2-3) 

We denote conditions at the unperturbed location by subscript 0, conditions at a perturbed 
location by subscript 1, and distinguish conditions inside the perturbed ring from those 
on the outside by a prime. Now concentrate on a portion of the loop that moves from ro 
to ri = r + Ar. Since the change is slow, hydrostatic equilibrium in the z— direction is 
maintained locally as a function of r, implying 



d_ 

dr 



I \2 
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8tt 



d_ 
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+ Pi 
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(2.4) 



where the final equality comes from eq. (2.2). Moreover, if we integrate eq. (2.4) over 
radius the constant of integration must vanish so that 

<s£ (2 ,) 

5 



In general, hydrostatic equilibrium cannot be satisfied simultaneously in both the 
r— and z— directions. There will either be a restoring force in the r— direction, which 
moves the loop back toward its original location, or a destabilizing force which augments 
the perturbation. When fc ^ 0, one must include the magnetic tension force due to 
the perturbed B z . For m = 1 perturbations that are periodic in the z— direction with 
wavenumber k, the force due to the perturbed B z can be written as the sum of two terms, 
— (2kbB/r + k 2 b 2 )Ar/4ix. The quadratic term (oc (kb) 2 ) results from the distortion of B z 
into B r , giving rise to a tension force in the r — z plane which always acts as a restoring 
force. The term linear in kb results from the distortion of B z into a combination of B r and 
B^. The resulting tension force (in the r — <j) plane) can either stabilize or destabilize the 
perturbation, depending on whether kb is positive or negative. If we visualize the m = 1 
mode as a helix wound about the z— axis, then kb > corresponds to a helix wound in 
the same sense as the magnetic field. The increased twist in B z then increases the tension 
in the r — (p plane, yielding a stabilizing effect. Conversely, a mode with kb < is wound 
counter to the magnetic field. The twist in B z then decreases the net tension, helping to 
destabilize the mode. 

The condition for stability is 



d_ 

dr 



t\2 



(gj) 
8tt 



(B'V Ar 1 

+ + (2kbB + k 2 rb 2 ) -=— \ Ar > 0. (2.6) 

Anr 47rr J 



We use r in the term for the perturbed magnetic hoop stress because the radius of the 
loop is unchanged. Eliminating the partial derivative through the use of eq. (2.4) and 
noting that 

(B[f Bl _(B* + 2 £A Ar 



47rro 4-7rri \4-7rr 2 r dr 
where we have used eqs. (2.5) and (2.3) (and have dropped subscripts as permitted, to 
first order), we obtain 



B 2 8n 

— + + (2kbB + k 2 rb 2 ) 

r or 



(Ar) 2 > (2.8) 



as the condition for stability . Substituting for the pressure gradient from eq. (2.2) and 
noting that (Ar) 2 is positive-definite, we finally obtain the necessary condition for stability 
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of m = 1 modes: 

d^B 1 , b 1 / b\ 2 

-d^ < -2 +kr B + 2[ kr B) • (2 ' 9) 
Thus, any region of a Z-pinch in which B increases with radius or decreases more slowly 
than r -1 / 2 should be subject to m = 1 instability for 1 <tikr <ti B /b. We verify this result 
using the full dispersion relation in § 3. 

Next consider m = modes, which correspond to changing the radius of a flux loop 
(from ro to r\ = tq+ Ar) while keeping its center fixed. Pressure balance in the z— direction 
again requires that eqs. (2.4) and (2.5) be satisfied. The stability condition is then 



d_ 

dr 



(B',) 2 .1 (SO 2 k 2 rb 2 Ar 
+ Pi ' 



8n 



+ + ^ }Ar>0. (2.10) 

47rri Anr r " 



Note that we now use T\ in the hoop stress term (instead of ro) because the radius of the 
loop has changed, and that only the quadratic term in the B z tension force appears for 
m = perturbations. Using eqs. (2.2) and (2.4) and defining AS = B[ — B , we can 
manipulate eq. (2.10) into the form 

rflnS AB/B 1 / b \ 2 , n „. 

<xr + d^d • (2-11) 



<ilnr Ar/r 2 \ B 

To evaluate (AB / B)/(Ar /r), we define Ap = p[ — p and use eqs. (2.2) and (2.5) to 
show that 

4,Ap = -B 2 + ^) • (2-12) 

We must then consider the cases of finite and vanishing b separately. If b = identically, 
then flux loops are truly closed. The gas pressure responds to changes in the volume of 
the loop, V, according to an adiabatic law, 

Ap AV ( AA Ar\ 

where A is the cross-sectional area along the loop. From magnetic flux conservation, we 
also have 

— = ~- (2.14) 
B A K ' 
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Combining eqs. (2.12), (2.13), and (2.14) with eq. (2.11), and adopting the usual plasma 
'beta' parameter f3 = 8tcp/B 2 , we obtain the stability criterion 

d In B 70-2 

~dh^ < 7^+2' {2Ab) 

Surprisingly, the stability threshold with b ^ is quite different from eq. (2.15), even if 
b is vanishingly small. In this case, the flux loops are connected to form a cylindrical helix. 
Near marginal stability, fluctuations are extremely slow so there is time for the pressure 
to equilibrate along the helix. Therefore Ap = 0, eq. (2.12) gives (AB/B)/(Ar/r) = —1, 
and the stability criterion follows immediately from eq. (2.11). In § 3 we verify both cases 
of the m = stability threshold, and assess their physical implications in later sections. 

3. DERIVATION OF DISPERSION RELATION 

3.1. Linearized Equations 

The relativistic momentum equation for a magnetofluid can be written in the form 

r2 (^f)S-Vv) + V P+ j|-.E-i^ = (3,, 

(e.g., Weinberg 1972), where T 2 = (1 — w 2 /c 2 ) -1 , p and p are the proper pressure and total 
energy density, respectively, p e is the charge density, and the electromagnetic fields and 
current density are represented by their usual symbols. In this paper we are concerned with 
perturbations about time-independent, static equilibria, and therefore we treat the velocity 
v as a first-order quantity and T — 1 as second-order. By the usual arguments of MHD, the 
electric field E is first-order, and since p e = V • E/47T the electrostatic force term is second- 
order. Distinguishing zeroth from first-order quantities by the corresponding subscripts, 
we obtain the following first-order momentum equation: 



Po\ dw ji x B j x Bi 

— = -Vpi H H , 3.2 

at c c 



which must be solved subject to the perturbed Maxwell's equations: 



„ „ 1 <9Ei 47TJ! 

VxB 1 = -— -i + ^i 3.3 

c at c 

Vx El = -i^ (3.4) 
V-Bi = 0; (3.5) 



Ohm's law: 



_ v x B 

Ej + = 0; 



and an equation of state, which we take to be 

dp 



dt 



I = _( v • V)p - 7Po(V • v), 



(3.6) 



(3.7) 



where 7 is the adiabatic index. Note in equation (3.3) that we retain the displacement 
current <9Ei/<9£, which is usually dropped in nonrelativistic MHD analyses. 
Using eq. (3.6) to eliminate Ei from equations (3.3) and (3.4), we obtain 



which includes the displacement current, and 

0Bi 



(3.8) 



dt 



= V x (v x B ), 



(3.9) 



the usual flux-freezing equation. We now substitute for the current densities in eq. (3.2) 
using eq. (3.8) for ji and 47rj /c = V x B , obtaining (after some manipulation) 



j^Bo + j^Bi = 1 v)Bq + _ v(Bq B 

c c 4n 



(3.10) 



Next we substitute eq. (3.10) into eq. (3.2), differentiate the entire equation with respect 
to time, and substitute eq. (3.7) to eliminate p\. After some rearrangement of terms, the 
resulting equation can be written 



P0 + n + BlY^ 



c 4tvc 2 J dt 2 4itc 2 



Br 



a 2 v 

dt 2 



V[(v- V)p +7Po(V- v)] 



+ 



1 

47T 



0Bi 



V B + (B • V) 



9Bi 



- V B • 



dB 1 
~dt 



(3.11) 



Note that the displacement current gives rise to an effective inertia associated with the 
magnetic energy density, for velocity perturbations which are transverse to the unperturbed 
magnetic field. One can derive a homogeneous linear equation in v by substituting for 
dUi/dt from eq. (3.9), but before taking this step we will first adopt some simplifying 
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assumptions about the nature of the underlying equilibrium and of the perturbations of 
interest. 

3.2. Quasi-Local Perturbations About Toroidal-Field Equilibria 

As in § 2, we specialize to perturbations about cylindrically symmetric equilibria (in 
coordinates {r,z,4>}) satisfying 

Po =Po(r), B = B Q (r)4> + bz, (3.12) 

where b is a constant. Such equilibria must satisfy eq. (2.2). We further specialize to short- 
wavelength (quasi-local) modes for which all perturbed quantities can be approximated to 
have the form 

^i(m(f>+kz-\-lr — ult) -^g^ 

The main requirement for such modes (assuming d\nB/d\nr ~ 0(1)) is that fc 2 r 2 , l 2 r 2 3> 
1. Since the background equilibrium depends on radius, true normal modes of the system 
have the more general form f ( r ') e t ( m 4>+kz-ujt) ^ rp Q self-consistency of our approx- 

imation scheme, we derived the exact global dispersion relation for the case 6 = 0, and 
showed explicitly that the exact equations reduce to our approximate equations under the 
assumptions adopted below. 

The quasi-local equations admit two classes of modes: "high-frequency" modes, the 
relativistic versions of magnetosonic and Alfven waves, and "low-frequency" modes, which 
include the instabilities we wish to study. To select the latter, and to determine which 
terms to retain in the quasi-local approximation, we adopt to following orderings, the 
self-consistency of which we have verified a posteriori: 

A;V>l + m 2 , 0(k)~0(l), (3.14) 
|^ 2 |<A; 2 c 2 , kv z + lv r ~ O (y) , (3.15) 

6<S , bl,bk~o(J^J. (3.16) 

To simplify notation, we will henceforth drop the subscript on unperturbed quantities. We 
then define the relativistic enthalpy, 

e = pc 2 +p, (3.17) 
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and the following dimensionless parameters: 



v 



e u 2 r 2 



7P & 



r dB b 8ivp 

a =B^ ^ kr B' P =W 



(3.18) 



Using the approximations and parameters defined above, the components of dTSi/dt 
in eq. (3.9) are given to lowest order by 



<9B] 

~dt 



k fdB 1 \ B . 



(3.19) 



(<9B \ B 
-^j = -[(l-aK-r(V-v) + z(m + r/)^], (3.20) 

where we have used eq. (3.15) to justify setting v z = ~(l/k)v r to lowest order in eq. (3.19), 
and have substituted kv z + lv r = — z(V ■ v) — mv^/r + iv r /r in eq. (3.20). 

We now manipulate the momentum equation to derive the dispersion relation. Writing 
eq. (3.11) in the form 

a = V/i + F, (3.21) 

where a represents the acceleration terms and h contains the gas pressure forces, we find 
the following expression for the components of the magnetic force, F, again keeping terms 
to lowest order: 



F r = 



B 

4%r 



I 2 



B r - (a + 1)B 4 



d(rB q 
dr 



F z = —ikr 
B 



BBtf, 
4/rr 



Anr 



(a + l)B r + i7]B^ , 



(3.22) 

(3.23) 
(3.24) 



where the dotted quantities represent the time derivatives (3.19) and (3.20). We will 
ultimately obtain the dispersion relation by evaluating the combination 



da z da r 2 / p 

~S IT =lkw \P + — + 

or oz \ c 

dF z dF r 



B 2 \ f P_ 

2 ' 4ttc 2 { + k 2 



dr dz 

. — I 2S - i(m + rj) ( 1 + L 1 B, 



4nr 



k 2 



(3.25) 
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where we have used the ordering adopted above to retain terms only to lowest nonvanishing 
order. Note that the gas pressure terms cancel in this combination. Dividing through by 
ikB 2 1 (4irr 2 ) , substituting for B r from eq. (3.19), and rearranging terms, we obtain 



rB 4 
B 



(3.26) 



where v is the dimensionless frequency defined in eq. (3.18). 

We use the z— and 0— components of eq. (3.21) to eliminate B$ (or, implicitly, 
V • v and v$) in favor of v r . In the z— component, we may neglect the acceleration terms 
compared to the gas pressure and magnetic forces (just as, in our heuristic treatment 
[§ 2], we assumed magnetohydrostatic equilibrium in the vertical direction). Using the 
equilibrium condition written in the form 



dr Anr 
to eliminate the pressure gradient, we find (after some algebra) 

( 7 /3\ 

i(m + 77)^0 = I 1 + — I r(V ■ v) - 2v r 
This allows the right-hand side of eq. (3.26) to be simplified to 

^ = |r(V.v)-(« + lK 



(3.27) 



(3.28) 



(3.29) 



Using eqs. (3.27) and (3.28), the </>— component of the momentum equation can be simplified 
to 

— Wff) = i(m + rj)r(V ■ v). (3.30) 
Using eq. (3.28) a second time to eliminate v^, we obtain 



(m + rj) 2 -v 1 + 



7/9 



r(V • v) = — 2vv r 



(3.31) 



Finally, eliminating r(V • v) in eq. (3.29), we are able to cancel factors of v r in eq. (3.26) 
to obtain the dispersion relation: 



/ |2\ (m + if) 2 - v [\ - 
V k2 J (m + 77) 2 -v(l + 3f) 



+ a = 0. 



(3.32) 
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Equation (3.32) is the principal result of this section. 

Note that eq. (3.32) can be written as a quadratic in u, f{y) = Av 2 + Bv + C = 0. It 
is easy to show that f{v) must change sign between its limit at v — > ±oo and its value at 
v = (m + r]) 2 /(l + 7/3/2). This implies that B 2 - 4AC > 0, i.e., that f(v) always has two 
real roots. Since v is effectively a; 2 , the modes are always either purely growing or purely 
damped as, of course, they must be for a nondissipative hydromagnetic system. 

4. INTERPRETATION OF DISPERSION RELATION 

4.1. Instability Threshold 

Unstable modes have imaginary uj, implying v < 0. We therefore seek values of 
l 2 /k 2 > 0, such that f{y) = for v < 0. Note that for l 2 /k 2 -> oo, eq. (3.32) implies that 
v — > (m + i]) 2 /(jp/e + 7/3/2) > 0, so the modes are stable in this limit. Unstable modes 
exist only for a finite range of l 2 /k 2 : 

0<^<2 . a + \- -l, (4.1) 

where the latter portion of the inequality comes from setting v = in eq. (3.32) with 
(m + rj) taken to be finite. The necessary and sufficient condition for the existence of 
instability is then 

dlnB 1, , . , nS 

which agrees exactly with the results of our heuristic analysis for both m = 1 (eq. [2.9]) 
and m = (eq. [2.11] with (AB/B)/(Ar/r) = -1). If m + rj = identically, the limit 
of eq. (3.32) as v — > yields eq. (2.15), also in agreement with our heuristic analysis for 
m = b = 0. 

4.2. Fastest Growing Modes 

The fastest growing modes are those with the most negative values of v. Recalling 
that v = (e/7p)a; 2 r 2 /c 2 , one sees from eq. (3.32) that growth rates depend on the radial 
wavenumber, /, only through the squared ratio of radial to vertical wavenumber, 5 = l 2 /k 2 . 
Growth rates also depend on k through r\ = kr(b/B). 
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As one might surmise from eq. (4.1), the growth rates of unstable modes increase with 
decreasing /, a result which can be verified by differentiating eq. (3.32) with respect to 5, 
holding k and m fixed. Defining a = (m + r]) 2 , d = 7/3/2, and e = 'jp/s, we find 



1/ r x 2ad 

(e + d)(l + 6) + 



2 V ;v ; [a-(l + d)z/] 2 



^ = i[a-(e + d)i/]. (4.3) 



The quantity in square brackets on the left-hand side is positive-definite for all v, while 
the right-hand side is positive-definite for v < 0. Therefore, du/dd > for i/ < 0. This 
means that the fastest growing modes occur for 5 <C 1, i.e., for vertical wavelengths much 
shorter than radial wavelengths. We will henceforth set l 2 /k 2 « in eq. (3.32), which now 
takes the form 

* [(e + cQ ,- a] + _|^_ + a + 1= o, (44 » 

Growth rates depend on and the mode number, m, through the combination a = 
(m + rj) 2 . Although m must be an integer or zero, a can take on any positive value, and 
unstable modes exist for < a < 2(a + 1). Differentiating eq. (4.4) with respect to a, we 
obtain 



1 . 2ad 
(e + d) + 



dv 1 2du , . 

!hi = 2 + [a- (l + d)v] 2 ' ( ' ' 



_2 V ' [a-(l + d)v] 
The quantity in square brackets on the left-hand side is positive-definite, but for v < the 
right-hand side may be either positive or negative. As a approaches 2(a + 1), v approaches 
zero and dv/da is positive. This ensures that the growth rate increases with decreasing o 
near the stability threshold. As a decreases toward zero, two behaviors are possible: 1) v 
is monotonic in a with the largest growth rate at a — * 0, or 2) v exhibits a local minimum 
at finite a. To see that there can be at most one local minimum and no local maxima, 
differentiate eq. (4.5) with respect to a and set dv/da = 0, which shows that d 2 vjda 2 > 
for v < 0. 

To determine which of the two behaviors applies, evaluate eqs. (4.4) and (4.5) at a = 
and eliminate v between them. If dvjda < at a = 0, then the maximum growth rate 
occurs at finite a. This condition is equivalent to an upper limit on a, 

d-1 2d(e + d) tl , 
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The inequality (4.6) is satisfied for a wide range of plausible equilibria, where: 1) p decreases 
monotonically outward, and 2) the toroidal flux surfaces do not contain a neutral sheet (a 
surface where changes sign). In such equilibria, a > — 1 everywhere, and a — > — 1 for 
d — > (low-/3 limit). If we approximate this class of equilibria by the analytic model 



a 



aod — 1 
d+1 



(4.7) 



then eq. (4.6) is satisfied for ao < 3. Note that eq. (4.7) with ao = 1 is satisfied identically 
for purely toroidal MHD equilibria established by adiabatic flow with = 0. These 
equilibria form a two-parameter family that has been studied by Begelman & Li (1992) as 
a model for the evolution of the Crab Nebula, and may also be applicable to jets as they 
propagate and expand laterally. Note also that models with ao < 1 have divergent current 
density on the axis. 

To find the maximum growth rate for systems satisfying eq. (4.6), we set the right- 
hand side of eq. (4.5) equal to zero, and solve the resulting equation simultaneously with 
eq. (4.4). The minimum value of v, 

(a + l)(l-ey 1/2 ' 



8d 



1 - 



(q+l)(l-e) 
4d 



1 - 



2d 



is attained at 



a r 



a + l + -(e + 2d+ l)v r , 



(4- 



(4.9) 



We chose the solution of the quadratic leading to eq. (4.8) that gives positive a m . 

4.3. Useful Approximations and Qualitative Summary 

The dispersion relation (4.4) is particularly simple to analyze in the limit of a weak 
magnetic field, (3 » 1 (d> 1). To lowest order, the result is 



v = d 1 1 + a - a - y/(a - l) 2 + 4a 

Evaluating equations (4.8) and (4.9) in the weak-field limit, we find 

(a + 1) 2 

Ad 



(4.10) 



(4.11) 
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with maximum instability occurring at 



(4.12) 

Equations (4.10)-(4.12) give remarkably accurate results over a much wider range of insta- 
bility than the asymptotic regime where the approximations are rigorously justifiable, i.e., 
where the field is weak. Because of the structure of the full dispersion relation, eq. (4.10) 
automatically gives the correct condition for v to vanish, which occurs at rather small 
values of f3 (e.g., at (3 = 2/(37) f° r a = a o = !)• Thus, the approximation is "pinned 
down" outside its nominal range of validity, a feature which could account for its useful- 
ness. The accuracy of the approximation increases with increasing a and decreasing ct . 
At a = ao = 1, it is accurate to within about 15% for all d (and gets better for higher a 
and lower ao), but it is rather poor near the stability threshold, for a — > 0. 

Evaluation of the dispersion relation for different parameter combinations reveals other 
important trends, m = 1 modes (with \rj\ <C 1) are more unstable than m = modes in 
the weak-field limit, provided that a < 5/2. In regions where the field is relatively stronger 
(P ^ 1), m = modes are relatively more important, and may begin to dominate for j3 
below some threshold value /3 t h > 2/(37). For ao close to 1, the crossover in growth rates 
occurs close to the stability threshold for the m = 1 mode. As we have noted, in the 
special (but physically motivated) case a = 1, there is no m = instability at all. Note, 
however, that the a = 1 modes are stable for a < — 1/2, whereas instability of a = modes 
can persist until a — > — 1. 

Despite the relationship of m — 1 and m — growth rates, there is no trend toward 
increasing instability with higher values of a (or m). Indeed, eq. (4.12) implies that the 
fastest growing modes in the weak-field limit always occur for a < 1. 

Finally, we note that the parameter e, which ranges between (c s /c) 2 <^ 1 in the 
nonrelativistic limit (where c s is the speed of adiabatic sound waves) and 1/3 in the 
ultrarelativistic limit, has relatively little effect on the growth rate for plausible values 
of d, a, or ao- 

To summarize, the toroidal-field equilibria under consideration are most unstable to 
m = 1 and m = modes, m = 1 modes dominate in the weak-field limit when ao < 5/2; 
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a r 



(a + 1) 



1 - 



(a + 1) 



when a is close to 1, they dominate throughout. (Of course, m > 1 modes can behave 
like lower-m modes if r] is negative with a large absolute value, but we will not discuss this 
possibility here.) For our purposes, eq. (4.10) or (4.11) will suffice to give rough estimates 
of the growth rate for the most unstable modes. Setting a = 1 and restoring dimensions, 
we have 

r 2 B 2 

u 2 „ — . (4.13) 
The corresponding growth timescale is 

t 3 = |^|- 1 ~(2/3) 1/2 - (4.14) 



c 



in the ultrarelativistic limit (s = 4p) , and 

t, =~ (4.15) 

where va is the Alfven speed, in the nonrelativistic limit. 

4.3. Nonlinear Development 

The growth timescales derived above are sufficiently short, compared to the time 
available in systems like the Crab Nebula and jets, that the nonlinear development of the 
instabilities must be considered. Given the problem's complexity, a qualitative analysis 
will have to suffice. Even at this level, it is clear that the m = and m = 1 modes will 
produce dramatically different effects. 

Consider the m = modes first, supposing that these grow faster than the m = 1 
modes. Their initial effects will be to create gradients of in the z— direction, without 
destroying the concentric field topology. Equilibrium generally cannot be established under 
such conditions, suggesting that reorganization of the field's radial structure will continue 
well into the nonlinear regime. What is the likely outcome of this reorganization? 

If the flow is adiabatic (even if the adiabats are different on different streamlines), then 
the system should evolve toward the family of equilibria derived by Begelman & Li (1992). 
These equilibria satisfy eq. (4.7) with ct = 1 5 an d are neutrally stable to m = r\ = modes 
at all radii. We therefore conjecture that m = modes, where dominant, will rearrange 
the field in such a way as to approach marginal stability. Once this happens, we need only 
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consider the m = 1 modes, which will be unstable for all (3 > 2/(37) (= 1/2 (2/5) for an 
ultrarelativistic (nonrelativistic) equation of state). (If B z 7^ 0, then m = modes will 
still be slightly unstable, with growth rates ~ 0(rj); however, the phase space of rapidly 
growing modes will be greatly reduced compared to the m = 1 modes, and the latter 
should dominate.) 

Our heuristic derivation of m = 1 instability in § 2 seems also consistent with contin- 
ued growth into the nonlinear regime. Recall that modes with larger radial wavelengths 
(smaller radial wavenumber /) are more unstable than modes with shorter radial wave- 
lengths, for fixed vertical wavenumber k. Thus, relatively large sideways excursions of 
narrow field loops are possible before strong nonlinearity sets in. The geometry of the 
unstable motion is reminiscent of a stack of coins coming apart. 

Once the displacement of a loop of field approaches and exceeds its initial radius, the 
loop no longer encloses the axis of symmetry and the dynamics of the instability becomes 
qualitatively different. Conditions would then seem to favor a significant reduction in mag- 
netic field energy through contraction and reconnection of flux loops. Free of the inhibiting 
effects of a concentric field structure, closed flux loops would tend to contract, adiabatically 
converting magnetic field energy into particle energy. Neighboring flux loops might also 
reconnect, turning magnetic into particle energy non-adiabatically. There seems to be no 
reason to expect a concentric field structure to be restored, or its destruction inhibited, by 
nonlinear effects. Modeling these effects is beyond the scope of this investigation. However, 
we make the following conjectures based on the nature of the unstable modes and their 
probable development: 

1) Regions of concentric toroidal field with (3 > 2/(37) are subject to m = 1 instabil- 
ity. This will include the axis and regions nearby, unless the longitudinal (B z ) field 
component becomes dynamically significant in these zones. 

2) Concentric field structure in the region affected by instability is not preserved. Once 
m = 1 instability becomes well-established, it will probably wash out significant mean 
field gradients in the affected region. Depending on the efficiency of reconnection, the 
reorganized field may retain considerable small-scale structure, with significant local 
gradients. In particular, if the field structure evolves through axisymmetric expansion 
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of a medium with an initially weak field — a plausible scenario in jets and in the 
Crab Nebula — then the field gradient may never become strong enough (f3 > 2/(37), 
d log B/d log r < —1/2) to shut off the instability. The reorganized field might remain 
predominantly aligned with the toroidal plane, making it appear (from polarization 
measurements) that the field remains well-organized. 

Since the Z-pinch is unstable down to very short wavelengths (and the presence of a 
small B z may even enhance instability depending on the helicity of the mode), inclusion 
of a finite resistivity may lead to some dissipation of magnetic energy and nonadiabatic 
behavior even in the linear regime, as well as to damping of the instability on the smallest 
scales. 

5. APPLICATIONS 

As argued in the Introduction, the magnetic field structures of astrophysical jets and 
the interiors of plerions tend to become dominated by concentric loops of toroidal field. 
Even if magnetic stresses dominate over pressure forces at large cylindrical radii, there must 
be a range of radii within which the condition (3 > 2/(37) holds. This region is susceptible 
to the m = 1 instability discussed in this paper. As loops of field in the unstable zone are 
dislodged from concentric alignment by the instability, field loops from outer, stable zones 
may shrink to take their place. But these too will then be subject to the m = 1 instability. 
Thus we claim that instability will substantially destroy the organized field structures of 
astrophysical Z-pinches, provided that the unstable modes have time to grow into the 
nonlinear regime. Below we estimate the growth criterion for generic conditions in jets, 
and for the best-studied plerion, the Crab Nebula. In each case, operation of the instability 
has important consequences for our understanding of physical conditions in these objects. 

5.1. Jets 

The growth times evaluated in § 4 are measured in the frame of the material moving 
down the jet. Thus, if we examine a region of a jet at distance R from the source, the 
instability will grow provided that Tt g < R/v, where v is the jet flow speed and the bulk 
Lorentz factor, V = (1 — f 2 /c 2 ) -1 / 2 , corrects for time dilation. If r is a characteristic 
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transverse scale across the jet (i.e., the jet's width), then the condition for growth of the 
instability in a nonrelativistic, supersonic jet is 

-<^ = — , (5-1) 
R v Ma 

where Ma is the Alfven Mach number. In a relativistic jet, it is convenient to express 
the growth criterion in terms of the ratio of Poynting flux to kinetic energy flux, the 
magnetization parameter (Michel 1969) 

ffS I? (8 ' 2) 

where B and e are measured in the fluid frame. The condition for growth of instability in 
a relativistic jet is then 

r a 1 / 2 

u < — < B ' 3 » 

Since most jets are observed to expand laterally as they propagate away from the source, 
the ratio r/R can be identified roughly with the rate of expansion of the jet. If the jet's 
interior is in causal contact with its surroundings — e.g., if it is pressure-confined — then 
the rate of expansion cannot exceed the Mach angle, i.e., 

" <U-Y /2C - (5-4) 



R r \e 

for a region dominated by gas pressure, where we have neglected numerical factors of order 
unity to emphasize the basic scaling. Observations of many jets indicate interactions with 
their environments, implying that eq. (5.4) is probably satisfied (e.g., Bicknell & Begelman 
1996). 

Equation (5.4) alone is not sufficient to guarantee that the instability growth condition 
(eq. [5.1] or [5.3]) is satisfied. The upper limit on r/R exceeds the instability growth limit 
by a factor ~ /3 1 / 2 , which is larger than one for a gas pressure-dominated jet. Yet there 
are two conditions under which a jet flow will evolve toward satisfying the instability 
growth criterion with increasing distance from the source, even if it starts out stable 
under the criterion. First, a pressure-confined jet, propagating down a pressure gradient 
shallower than r -2 , will become increasingly collimated (Begelman, Blandford, & Rees 
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1984). The ratio r/R decreases with R, while the right-hand side of the growth criterion 
remains constant (for matter-dominated jets) or increases (for relativistic flows with p 3> 
pc 2 ). Second, the growth criterion will be satisfied automatically if the mean (5 in the jet 
approaches unity. This can occur even at constant r/R, since the gas pressure in adiabatic 
expansion decreases more steeply with r than does B 2 . 

The conditions under which jets form are not well enough understood to make detailed 
predictions about the radius at which the instability will develop. We note that its onset 
may take place far from the jet's source and could lead to the dissipation of a large 
fraction of the magnetic energy into energetic particles. There is observational evidence 
for the onset of dissipation at some distance from the source of AGN jets, in the form of 
synchrotron "gaps" near the jets' bases (Bridle & Perley 1984). An analysis of synchrotron 
emission along the jet in M87 also indicates that much of the particle acceleration takes 
place at 1-10 parsecs (10 3 — 10 4 Schwarzschild radii) from the central black hole (Heinz & 
Begelman 1997), i.e., in a region quite distinct from the jet formation zone at <^ 0.01 pc 
(Junor & Biretta 1995). 

5.2. The Crab Nebula 

To a first-order approximation, the Crab Nebula can be described as a bubble inflated 
by a mixture of relativistic particles and toroidal magnetic field injected by the central 
pulsar in the form of a relativistic wind (Rees & Gunn 1974; Kennel & Coroniti 1984a; 
Begelman and Li 1992). A standing reverse shock is inferred to lie about 10% of the 
way from the pulsar to the edge of the synchrotron nebula, at the outer boundary of the 
"underluminous" region surrounding the pulsar (Kennel & Coroniti 1984a, and references 
therein) and just interior to the time-variable optical "wisps" (Scargle 1969; Hester et 
al. 1995). 

Outside the shock, hydrostatic equilibrium is maintained by a balance between rel- 
ativistic particle pressure and the combination of magnetic pressure and tension. The 
shocked wind must decelerate from J> c/3 just outside the shock to the nebular expansion 
speed of ~ 1500 km s _1 at the bubble's edge. Flux-freezing and gasdynamical arguments 
then lead to the following picture for the magnetic field structure. In the inner regions 
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the field strength increases linearly with radius, while the dominant particle pressure is 
approximately uniform. If the bubble extends beyond the radius at which the magnetic 
pressure becomes comparable to the particle pressure, then the structure would increas- 
ingly resemble that of a pure magnetic pinch, with B declining as r _1 and gas pressure 
becoming relatively unimportant. 

The inference that the magnetic field strength is an increasing function of radius in the 
inner nebula has been used to constrain the form in which energy is injected by the pulsar 
wind. If the ratio of magnetic pressure to gas pressure were larger than a few percent just 
outside the reverse shock, then the outer nebula would be strongly pinched and therefore 
highly elongated, in contradiction to observations (Rees & Gunn 1974; Begelman & Li 
1992). Moreover, the reverse shock would then be pushed inward from its observationally 
inferred radius, perhaps collapsing altogether (Rees & Gunn 1974; Emmering & Chevalier 
1987). Thus, in order for the Rees-Gunn model to satisfy boundary conditions at both 
the reverse shock and the nebular boundary, the magnetic pressure must not become 
significantly larger than the gas pressure anywhere in the nebula. Taking into account the 
amplification of the magnetic field due to compression at the shock front, this implies that 
the ratio of Poynting flux to kinetic energy flux (a) just inside the shock is no more than 
a few times 10 -3 . 

Begelman & Li (1992) confirmed this qualitative behavior in an exact magnetohydro- 
dynamic (MHD) model, and pointed out that the structure depends on cylindrical, not 
spherical, radius. Although local interactions between the relativistic fluid and the denser 
filaments of ejecta undoubtedly occur, and may partially be responsible for confining the 
nebula (Hester et al. 1996), the Rees-Gunn model seems to provide a robust description 
of the nebula's gross features. Thus one can view the Crab Nebula as a confined, cylindri- 
cal structure with both magnetic and particle pressure varying with cylindrical radius. In 
other words, to an excellent approximation the interior of the Crab Nebula should resemble 
a Z-pinch. 

Indeed, the Rees-Gunn model depends crucially on the assumption that the magnetic 
field retains a concentric toroidal geometry as the nebula expands. Is this realistic, given 
the instabilities to which Z-pinches are susceptible? To address this question, consider 
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the unique, two-parameter family of cylindrically symmetric MHD equilibria found by 
Begelman & Li (1992). These equilibria are characterized by the balance between par- 
ticle pressure and the forces due to a toroidal magnetic field, and subject to the dual 
constraints of magnetic flux- freezing and adiabatic flow along streamlines. Although the 
specific entropy of the gas is allowed to vary arbitrarily from streamline to streamline, the 
equilibrium structure turns out to be completely independent of the flow pattern. These 
models are ideally suited to describing the internal structure of the Crab Nebula and other 
plerions, at radii larger than the pulsar wind shock and in the limit that radiative losses 
have negligible effect on the equation of state. The models satisfy eq. (4.7) exactly, with 
a = 1, i.e., 

rflnr 7/3 + 2' 1 ' ' 

These equilibria are neutrally stable to m = modes, but are unstable to m = 1 modes for 
P > 1/2, where we have set the adiabatic index 7 = 4/3 to describe the ultrarelativistic 
fluid that is presumably injected by the shocked pulsar wind. 

At the radius of the pulsar wind shock, r s , the Rankine-Hugoniot conditions for a 
transverse MHD shock (Kennel & Coroniti 1984) allow us to deduce the downstream 
j3— parameter in terms of the magnetization parameter in the wind immediately upstream 
of the shock: 

A = 5^. (5.6) 

If we assume that f3 s is ^ 1 (but not necessarily as large as the values deduced by Rees 
& Gunn and their successors) , then flow just outside the shock is gas pressure-dominated. 
The postshock flow speed is then v ~ (c/3)(r/r s ) -2 , while (3 ~ (3 s (r/r s )~ 2 . Instability sets 



in at r,- where 



(see eq. [4.14]), giving 



^ ~ m 1/2 n/c 

V 



rsj 

r s 



l + (2/3 s ) 1 /2j V3 . (5.7) 



For the largest values of f3 s 150) derived using the Rees-Gunn model, eq. (5.7) gives 
ri/r s <> 3. 
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Thus the interior of the Crab Nebula should be highly susceptible to the local m = 1 
instabilities described earlier. If the instability reaches nonlinear amplitude as close to the 
pulsar wind shock as suggested by eq. (5.7), then very little of the nebula should possess 
the concentric field structure inferred by Rees & Gunn (1974). Loops of toroidal field 
deposited by the pulsar wind at r s will not simply expand about a common axis, but will 
be dislodged and "float away" before they expand very much at all. Instead of a field 
geometry dominated by nested loops of field a parsec or more across (i.e., comparable to 
the size of the radio synchrotron nebula, -R ne b) 5 we expect a chaotic structure with "loops" 
or field reversals common on scales ~ r s or less throughout the nebula (albeit possibly with 
a preferred field orientation — see below). Reconnection may be widespread, although we 
are unable to estimate its importance in dissipating magnetic energy. 

Destruction of the concentric field structure by m = 1 instability is not inconsistent 
with the high degree of polarization and preferred polarization direction observed in the 
Crab Nebula (Schmidt, Angel, & Beaver 1979; Michel et al. 1991). Recall that the most 
unstable modes have radial wavelengths much longer than their wavelengths along the 
cylindrical axis (§ 4.2). Since unstable motions occur mainly in the r — plane (at least in 
the linear regime), it seems plausible that B z remains small, even as and B r become 
randomized. If the rotation axis of the pulsar lies not too far from the plane of the sky, 
the randomized field would still produce a high degree of synchrotron polarization, with 
a direction preferentially parallel to the projected rotation axis. In fact, high resolution 
optical and radio observations of the nebula (Hickson & van den Bergh 1990; Bietenholz 
& Kronberg 1991), the latter corrected for foreground Faraday rotation, reveal a cellular 
structure, with a preferred direction overall but with fluctuations having a characteristic 
scale size of about 10" (= 0.1 pc at the assumed distance of 2 kpc). This scale size is very 
similar to the inferred value of r s , and supports our conjecture that the loops of toroidal 
field never grow beyond this scale before instability sets in. 

If our picture of a chaotic magnetic structure in the Crab is correct, then we must 
discard the dynamical argument used by Rees & Gunn (1974) and their successors to 
estimate a. To see this, it is useful to recall the "winding number" analogy introduced in 
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Rees & Gunn (see their eq. [4]). The magnetic energy inside the nebula satisfies 

2 

(5.8) 

where M is the number of turns made by the pulsar and £ is the mean size of a field loop. 
In the Rees-Gunn model, £ ~ -R n eb, since all of the loops expand concentrically. Rees 
& Gunn (1974) therefore implicitly set this factor equal to one, and did not include it 
explicitly in their formula. If the instability operates as we conjecture, then £ ~ r s and 
(£/R neh ) 2 ~ 0.01. To put this another way, the magnetic pressure would not increase with 
radius as r 2 . Therefore, to match the pressure and expansion rate of the nebula to its 
observed values, we can tolerate a value of a a hundred times larger than that inferred by 
Rees & Gunn (1974). 

6. CONCLUSIONS 

I have analyzed the local stability of axisymmetric perfect MHD equilibria in which 
gas pressure balances the pressure and tension of a toroidal magnetic field. This type of 
equilibrium — the "Z-pinch" of plasma physics literature — is germane to plerions, like 
the Crab Nebula, and probably applies to many astrophysical jets as well. The essential 
results of the stability analysis are not new — they can be found in any number of elemen- 
tary plasma physics or MHD textbooks, such as Bateman (1978) or Krall & Trivelpiece 
(1973). I present a complete derivation and discussion of the instability here because 1) 
the astrophysical applications require a relativistic treatment of the instability, which I 
have not found in the literature, and 2) I want to make the physics of these important 
instabilities more accessible to astrophysicists. 

Provided that they do indeed possess the assumed field geometry, jets and plerions 
are highly susceptible to internal kink instabilities, which probably destroy the concentric 
field structure. Unless jets possess a relatively coherent "backbone" of longitudinal (B z ) 
field, without field reversals on scales much smaller than the jet width, they are unlikely 
to develop any significant degree of magnetic self-collimation. However, the slow growth 
rate of the instability in a weak magnetic field implies that a well-ordered toroidal field 
structure could persist for a significant distance along the jet, depending on how tightly 
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collimated the jet is. Dissipation associated with the eventual onset of instability could 
provide a source of energy for particle acceleration at otherwise unremarkable regions of a 
jet far from its source. 

In the Crab Nebula we expect the instability to develop very close to the pulsar wind 
shock, at distances of only ~ 0.1 pc from the pulsar. Disruption of the concentric field 
structure would imply that the mean field strength is not amplified much by the expansion 
of the flow between the wind shock and the outer edge of the nebula. The problem of 
confining the concentric toroidal field, first noted by Rees & Gunn (1974), goes away, and 
with it goes their deduction that a — the ratio of Poynting flux to kinetic energy flux - 
in the pulsar's wind must be 1. 

The low value of a derived from the Rees-Gunn model has puzzled theorists for years. 
Perfect MHD models of pulsar winds predict that a cannot become much smaller than one 
(Michel 1969; Kennel et al. 1983; Arons 1992), unless highly contrived boundary conditions 
are satisfied (Begelman & Li 1994; Tomimatsu 1994). Attempts to explain the decline 
in a through dissipative effects have led to some ingenious but controversial proposals 
(Coroniti 1990; Michel 1994; Melatos & Melrose 1996). While analyses of the nebula's 
synchrotron spectrum and emissivity (Kennel & Coroniti 1984b; Hoshino et al. 1992) have 
lent support to a low-a wind, it has always been the dynamical argument that has seemed 
most incontrovertible. The proposal we have made — that the Rees-Gunn field geometry 
is physically unrealistic — does not allow us to determine the value of a in the Crab pulsar 
wind. However, it does eliminate the principal dynamical argument for concluding that 
a is small. Spectral arguments nonwithstanding , it would seem worthwhile to reconsider 
whether a might be as large as ~ 0(1), after all. 

I am grateful to Martin Rees for his valuable insights during the three-year gestation 
period of this paper, and for encouraging me to finish it. Ellen Zweibel critiqued an early 
draft and supplied advice on the plasma physics literature. This work has been supported 
in part by National Science Foundation grants AST-9120599 and AST-9529175. 



26 



REFERENCES 

Arons, J. 1992, in The Magnetospheric Structure and Emission Mechanisms of Radio 
Pulsars, IAU Colloq. 128, ed. J. A. Gil, T. H. Hankins, & J. M. Rankin (Zielona 
Gora: Pedagogical Univ. Press), p. 56 

Bateman, G. 1978, MHD Instabilities (Cambridge: MIT Press) 

Begelman, M. C, Blandford, R. D., & Rees, M. J. 1984, Rev. Mod. Phys., 56, 255 

Begelman, M. C, & Li, Z. 1992, ApJ, 397, 187 

Begelman, M. C, & Li, Z. 1994, ApJ, 426, 269 

Benford, G. 1978, MNRAS, 183, 29 

Bicknell, G. V., & Begelman, M. C. 1996, ApJ, 467, 597 

Bicknell, G. V., & Henriksen, R. N. 1980, ApLett, 21, 29 

Bietenholz, M. F., & Kronberg, P. P. 1991, ApJ, 368, 231 

Bridle, A. H., & Perley, R. A. 1984, ARA&A, 22, 319 

Chan, K. L., & Henriksen, R. N. 1980, ApJ, 241, 534 

Chen, F. F. 1974, Introduction to Plasma Physics (New York: Plenum) 

Coroniti, F. V. 1990, ApJ, 349, 538 

Eichler, D. 1993, ApJ, 419, 111 

Emmering, R. T., & Chevalier, R. A. 1987, ApJ, 321, 334 

Freidberg, J. P. 1982, Rev. Mod. Phys., 54, 801 

Heinz, S., & Begelman, M. C. 1997, ApJ, in press 

Hester, J. J., et al. 1995, ApJ, 448, 240 

Hester, J. J., et al. 1996, ApJ, 456, 225 

Hickson, P., & van den Bergh, S. 1990, ApJ, 365, 224 

Jackson, J. D. 1975, Classical Electrodynamics (New York: John Wiley & Sons) 
Junor, W., & Biretta, J. A. 1995, AJ, 109, 500 

Hoshino, M., Arons, J., Gallant, Y. A., k Langdon, A. B. 1992, ApJ, 390, 454 

27 



Kadomtsev, B. B. 1966, in Reviews of Plasma Physics, 2, ed. M. A. Leontovich (New 
York: Consultants Bureau), p. 166 

Kennel, C. F., & Coroniti, F. V. 1984a, ApJ, 283, 694 

Kennel, C. F., & Coroniti, F. V. 1984b, ApJ, 283, 710 

Kennel, C. F., Fujimura, F. S., & Okamoto, I. 1983, J. Astrophys. Geophys. Fluid. Dyn., 
26, 147 

Konigl, A., & Choudhuri, A. R. 1985, ApJ, 289, 173 

Krall, N. A., & Trivelpiece, A. W. 1973, Principles of Plasma Physics (New York: McGraw- 
Hill) 

Kruskal, M., & Schwarzschild, M. 1954, Proc. Roy. Soc, A223, 348 
Melatos, A., & Melrose, D. B. 1996, MNRAS, 279, 1168 
Michel, F. C. 1969, ApJ, 158, 727 
Michel, F. C. 1994, ApJ, 431, 397 

Michel, F. C, Scowen, P. A., Dufour, R. J., & Hester, J. J. 1991, ApJ, 368, 463 
Rees, M. J., & Gunn, J. E. 1974, MNRAS, 167, 1 
Scargle, J. D. 1969, ApJ, 156, 401 

Schmidt, G. R., Angel, J. R. P., & Beaver, E. A. 1979, ApJ, 227, 106 
Tomimatsu, A. 1994, PASJ, 46, 123 

Weinberg, S. 1972, Gravitation and Cosmology: Principles and Applications of the General 
Theory of Relativity (New York: John Wiley & Sons) 



28 



